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ABSTRACT 

The feasibility of using a multi-wavelength Mie-Raman lidar based on a tripled Nd:YAG laser for 
profiling aerosol physical parameters in the planetary boundary layer (PBL) under varying conditions 
of relative humidity (RH) is studied. The lidar quantifies three aerosol backscattering and two 
extinction coefficients and from these optical data the particle parameters such as concentration, size 
and complex refractive index are retrieved through inversion with regularization. The column- 
integrated, lidar-derived parameters are compared with results from the AERONET sun photometer. 
The lidar and sun photometer agree well in the characterization of the fine mode parameters, however 
the lidar shows less sensitivity to coarse mode. The lidar results reveal a strong dependence of particle 
properties on RH. The height regions with enhanced RH are characterized by an increase of 
backscattering and extinction coefficient and a decrease in the Angstrom exponent coinciding with an 
increase in the particle size. We present data selection techniques useful for selecting cases that can 
support the calculation of hygroscopic growth parameters using lidar. Hygroscopic growth factors 
calculated using these techniques agree with expectations despite the lack of co-located radiosonde data 
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Despite this limitation, the results demonstrate the potential of multi-wavelength Raman lidar 
technique for study of aerosol humidification process. 


1. INTRODUCTION 

The radiation balance of the Earth is strongly influenced by atmospheric aerosols of natural and 
anthropogenic origin (D’Almeida et al. 1991; Pilinis et al. 1995). For accurate modeling of 
corresponding radiative forcing knowledge of the vertical distribution of particle macro- and 
microphysical parameters is needed. During recent years ground-based and airborne lidars have 
become important tools for profiling tropospheric aerosols using either single or multiple wavelengths 
(Kovalev and Eichinger 2004). For quantitative studies of the optical properties of tropospheric aerosol, 
Raman lidars have proven to be most useful (Ansmann et al. 1992; Ferrare et al. 1998a,b; Ansmann et 
al. 2000; Turner et al. 2002). This lidar type measures elastically backscattered light simultaneously 
with Raman backscatter from molecules (nitrogen or oxygen), thus allowing independent calculation of 
particle backscattering and extinction coefficients without the need of critical assumptions about 
atmospheric parameters. Still, to invert the measured aerosol scattering coefficients to the particle 
microphysical parameters when only a single output wavelength is used, physical models with 
numerous a priori assumptions are generally needed. 

On the other hand, the application of multi-wavelength (MW) Raman lidar can be used to 
quantify the main particle microphysical parameters with fewer a priori assumptions. Information on 
the particle extinction and backscatter coefficients at multiple wavelengths is produced by such lidars 
and, with this infonnation, microphysical properties of aerosols can be retrieved through mathematical 
inversion (Muller et al. 1999; Veselovskii et al. 2002). Theoretical and experimental studies perfonned 


2 



during the last decade have demonstrated that the MW lidar technique is able to provide comprehensive 
information about aerosol microphysical parameters (Muller et al. 2004, 2005; Veselovskii et al. 2004, 
2005) although simulation studies have shown a decreasing sensitivity as particle size increases. These 
results demonstrate that the key to the successful retrieval of aerosol parameters is the joint use of 
aerosol backscattering and extinction coefficients, which can be accomplished through the use of 
Raman or high spectral resolution methods (Shipley et al. 1983; Liu 2002). Another important finding 
is that a simplified multi-wavelength Raman lidar based on a frequency-tripled Nd:YAG laser can be 
used to quantify the particle size distribution (PSD) and complex refractive index (Muller et al. 2001; 
Veselovskii et al. 2004, 2005). 

To study this technique further the lidar-derived parameters should be compared with 
comparable results from other instruments. One of the recognized instruments for retrieval of column- 
integrated aerosol properties is the robotic sun photometer (SP) that is used in the network called 
AERONET (Holben et al. 1998, Dubovik et al. 2002). The comparison of lidar and sun photometer 
data has been performed in previous publications (Muller et al. 2004) as a method to assess the 
performance of the multi-wavelength lidar technique. We also compare lidar-derived results with sun 
photometer, however we emphasize the height variation under different atmospheric conditions 
revealed by lidar versus the column-integrated quantities measured by the sun photometer. 

One of the key factors influencing the scattering properties of aerosols is their hygroscopicity. 
As the relative humidity (RH) rises, uptake of water vapor by hygroscopic aerosols increases the 
particle size, which generally leads to an increase in light-scattering (Hanel 1976). Numerous 
experimental studies of the relationship between aerosol scattering and relative humidity have been 
performed in the laboratory with aerosols of known chemical composition (Tang and Munkelwitz 
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1993; Fuzzi et al. 2001; Kumar et al. 2003; Svenningsson et al. 2006; Kreidenweis et al. 2006). In-situ 
experiments in the natural environment have also been performed using nephelometers (Kotchenruther 
et al. 1998, 1999; Carrico 2003) and differential mobility analyzers (Gasparini et al. 2006a, b) in tandem 
with humidity control systems. Under the right conditions, a lidar may provide another tool for 
studying aerosol hygroscopicity. The advantage of using a lidar (as opposed to laboratory or most in- 
situ studies) is that it can provide relatively continuous altitude-resolved measurements of aerosol 
properties without perturbing the aerosol or its surroundings. The obvious disadvantage is that, 
because the aerosol sampled by the lidar is not controlled in any way (e.g., source, state, etc), lidar 
studies of aerosol hygroscopicity need to be limited to occasions when the same aerosol type exists in 
at least a portion of a profile that is characterized by widely changing RH so that observed differences 
in aerosol properties can be attributed primarily to water uptake as RH increases. 

The enhancement of lidar backscatter due to particle swelling was first considered by 
MacKinnon, 1969. Since then numerous studies have investigated the relationship between aerosol 
backscattering and relative humidity (Ferrare et al. 1998b; Pahlow et al. 2006; Wulfmeyer and 
Feingold 2000; Feingold and Morley 2003; Raut and Chazette 2007). The water vapor mixing ratio in 
these more recent studies was derived from either lidar measurements, Raman (Ferrare et al. 1998; 
Pahlow et al. 2006) and differential absorption (Wulfmeyer and Feingold 2000), or was calculated 
under a boundary layer cloud deck suggesting a well-mixed boundary layer (Feingold and Morley 
2003). If the boundary layer is well mixed, then any change in backscatter in the vertical is considered 
to be due to uptake of water vapor by the particles. Assuming the atmosphere is saturated at cloud base 
and assuming that the potential temperature is constant with height, the relative humidity (RH) profile 
can be calculated (Pahlow et al. 2006; Wulfmeyer and Feingold 2000; Feingold and Morley 2003). 
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Using single-wavelength backscatter lidar data combined with thennodynamic assumptions of the 
mixing state of the atmosphere allows for the determination of the hygroscopic growth factor /(RH) for 
relative humidities close to saturation. The growth factor is calculated in terms of backscattering at a 
given RH, relative to some lower RH re f. The application of the multi-wavelength Raman lidar 
technique has an advantage over single-wavelength measurements because it increases the amount of 
available information allowing derivation of the particle size distribution and other physical properties 
along with the water vapor mixing ratio. 

The possibility of using the multi-wavelength lidar method for the study of cloud condensation 
nuclei was considered by Feingold and Grund in 1994, but experimental implementation of this method 
is more recent. The increase of particle size and extinction coefficient near the PBL top derived from 
multi-wavelength lidar measurements during the ACE 2 (Aerosol Characterization Experiment) field 
campaign has been reported (Muller et al. 2002; Ansmann et al. 2002). As inversion routines and 
retrievals improve and additional aerosol products become available (e.g., PSDs, refractive index), the 
use of multi-wavelength lidar to investigate aerosol hygroscopicity may provide new insights into the 
aerosol properties relevant for climate forcing in the atmosphere. 

The goal of this paper is to derive vertical profiles of aerosol physical parameters in the PBL 
from multi-wavelength lidar measurements under varying relative humidity conditions and to compare 
the column-integrated aerosol parameters with results from the AERONET sun photometer. To 
demonstrate the potential of the multi-wavelength technique to study aerosol hygroscopic growth, data 
selection techniques are used to assess the mixing of aerosols in the PBL and hygroscopic growth 
factors are calculated for the best of these cases. 
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2. Lidar system description 

The experiments described here were performed with a multi-wavelength Mie-Raman lidar 
developed at NASA/GSFC. The lidar is based on a Continuum 9050 laser with a 50 Hz repetition rate. 
The output powers at A,=355, 532 and 1064 mn were 17.5, 7.5 and 14 W respectively. The 
backscattered light was collected by a Meade LPX200 40-cm aperture Schmidt-Cassegrain telescope 
operated at 0.35 mrad field-of-view and inclined so that the elevation angle was 18 degrees. The off- 
vertical measurements obtained through this approach improve the extraction of parameters in the PBL 
by decreasing the altitude to which the lidar overlap function influences the retrievals. The collected 
light passes through an optical fiber to a receiving module, which includes an off-axis parabolic mirror 
for collimation. The spectral components of the collimated optical signal are separated by dichroic 
beam splitters and interference filters. In the current configuration the system is capable of detecting 
three backscattered and two nitrogen Raman signals at Lr=387 and 607 mn. The receiving module also 
has the ability to measure the water vapor and liquid water. Raman signals at 407.5 and 402 nm, 
respectively (Whiteman and Melfi 1999), although no liquid water data are analyzed here. The 
detection in the 355-532 nm range is perfonned with Hamamatsu R1924 PMTs; for the 607 nm 
channel a Hamamatsu H7422P-40 module is used; and backscatter at 1064 nm is measured with an IR 
enhanced Si APD (Licel APD-3.0) operating in analog mode. The outputs of the detectors are recorded 
at 7.5 m range resolution with a Licel data acquisition system that incorporates both analog and 
photon-counting electronics. 

3. Raman Lidar calculation techniques 
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The combined use of Rayleigh-Mie and Raman backscatters allows independent calculation of 
the particle backscattering [1 and extinction a coefficients at 355 and 532 nm (Ansmann et al. 1998). 
Computation of Pio 64 is performed from the backscattered signal using the following form of the lidar 
equation: 


p^) = A 


-C&"(z)+ft"(z)]exp 
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with the boundary condition of Pl( z re f) - 0 , where Pfiz) is the backscattered power, A>_ is a range 

independent calibration parameter, P “‘ , PI , a 7 , a l are molecular (Raleigh) and aerosol backscattering 
and extinction coefficients at wavelength X. The temperature dependence of Raman scattering has been 
neglected in these calculations due to the relatively mild change in atmospheric temperature over the 
altitude ranges involved (Whiteman 2003a, 2003b). The extinction profile at 1064 nm is derived by 
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extrapolating Raman measurements at 532 nm through the relationship — 7 - = ~r , where the Angstrom 
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coefficient k is calculated from the — ratio obtained from the lidar data and assumed constant over 

a 532 

the wavelength range. The reference distance z re f is chosen above the PBL top where aerosol 
contribution to backscatter is negligible. A typical lidar signal at 1064 nm is presented in Fig.l. It 
shows that for distance z re f>14 km (for height above 4.2 km) the lidar signal is mainly due to molecular 
scattering. 

The sensitivity of the calculated backscatter and extinction coefficients to the Angstrom 
coefficient, k, is generally mild (Whiteman 2003a) when using the standard Raman lidar technique 
(Whiteman et al. 1992). However, these previous results refer to soundings acquired vertically. For 
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small sounding angles with respect to the horizon, the optical depth of the sounding path is increased 
significantly and so is the sensitivity to k and its range dependence k(z). It should be mentioned that the 
Angstrom coefficient affects extinction and backscattering in different ways. The uncertainty in the 
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where k R - the wavelength of the Raman scattered signal, is 


determined by the fully propagated uncertainties in k(z). The backscattering coefficient 
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reference point. In our computations we use range dependent Angstrom coefficient k(z) derived from 
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the extinction ratio and neglect the spectral variation of k. From the computations performed for 

^532 

different k we conclude that the extinction is not very sensitive to variations in k, with corresponding 
uncertainties generally being less than 4%. However, for backscattering the use of a range-dependent k 
is essential. For some cases studied here, the uncertainty in the retrieval of p 355 at a height h=500 m can 
increase by as much as 10% if the range variation in k is not considered. When calculating Pio 64 , the 
approximation that k is spectrally independent in the 532-1064 nm range may appear too coarse. 
However, the aerosol extinction at 1064 nm is much lower than at 532 nm, so the backscattering at 
1064 nm is not very sensitive to spectral or range variations of the Angstrom coefficient. Based on 
these considerations, we estimate the uncertainty of both the backscattering and extinction calculation 
inside PBL to be less than 10%. Typical vertical profiles of backscattering and extinction coefficients 
at different wavelengths measured on August 22 are shown in Fig. 2. 


4. Retrieval of particle microphysical parameters 

The lidar- derived backscattering and extinction coefficients are then used to calculate aerosol 
microphysical properties through inversion. The main features of our retrieval algorithm, which is 
based on inversion with regularization (Twomey 1977; Tikhonov and V.Y. Arsenin 1977), are 
presented elsewhere (Veselovskii et al. 2002; 2004). The algorithm makes no assumption about optical 
data input errors. The only constraints on the pennitted refractive index and particle size distribution 
are that the refractive index is considered to be wavelength independent and that the concentration of 
the particles with radii below some r m i n and above some r max is zero, where the values of these radii are 
found in the process of inversion. Furthermore, in the retrieval we don’t consider particles outside the 
radii interval 0.05 pm - 10 pm, i.e. these are the limiting values of r min and r max . 

The algorithm was initially designed for retrievals at a single height layer, but in data analysis 
usually the whole vertical profile of particle parameters is desired. In principle, the profile can be 
obtained from the retrievals at every layer, though such procedure is time consuming and does not take 
into consideration the likely interdependence of aerosol properties between two closely spaced layers. 
In the new version of the algorithm we have introduced modifications to overcome these limitations. 

At any height layer the measured optical data g, (backscattering or extinction) are related to the 
size distribution / O') through the Fredholm integral equation: 

r 

max 

J K i (m,l,r)f(r)dr=g i i=l,..,L (2) 

Y ■ 

min 

where K t (m,r ,A) are the kernel functions for refractive index m and wavelength A calculated with Mie 
theory. Particle size distribution is approximated by the superposition of base functions d,( r ) as: 
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f(r) = y £ j C j (z)B j (r) ! 

j = i 

where c /( z )are the weight coefficients. The base functions have a triangular shape on a logarithmic- 
equidistant grid across the chosen radius interval [r m i n , r max ] in which the inversion is performed. 

The integral equation (2) transforms to a linear one as follows: 

‘i 

g,=Yu A fr (4) 

j = i 

where the elements of weight matrix A,, are calculated as: 

'max 

4 = { K i (m,A,r)B j (r)dr. (5) 

''min 

At any of / height layers the expression (4) can be written in matrix- vector form as: 

aV = g ; /=1 ,...,N l (6) 

or it can be expressed in more general form for the whole height range: 

AC = G (7) 

In the new version of algorithm the equation (7) is solved for all layers simultaneously by means of 
inversion with regularization: 

C = (A r A + yHr‘A r G (8) 

where y is the regularization parameter. Matrices A, H and vectors C, G have the following structure 

(Kolgotin and Muller 2008): 
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As shown in our previous publications (Veselovskii et al. 2002, 2004), the averaging of solutions 
in the vicinity of the minimum of discrepancy helps to stabilize the inversion. In the previous version 
of the algorithm the choice of the averaging interval was made for each layer separately. In the new 
version used here, all layers are processed simultaneously. The procedure includes the following steps: 
For every k - th solution the discrepancy pk is calculated (Veselovskii et al. 2002, 2004). At each 
height layer the solutions are ordered in accordance with their discrepancy pk, from small 
discrepancy to largest discrepancy. 

At every height layer we calculate the discrepacy (up for the first /?; averaged solutions. 

The total discrepancy for the whole height profile is then calculated as 
n t 
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To estimate the number of solutions to be averaged it is convenient to consider the dependence of /5 U)l 
on the individual pk, i.e. p tot (pk). Thus at every layer we sum the solutions with discrepancy smaller 
than pk . The analysis of the function p tot (pp permits the determination of the optimal solution set for 

all layers simultaneously where in general the optimum pk yields a stable p U>1 (pk) such that the 

increase of pk does not lead to a significant change in the derived microphysical parameters. 

The typical uncertainties of the retrievals studied here were estimated from modeling and from 
the scatter of the individual solutions to be 20% for volume density and effective radius, 25% for mean 
radius and 40% for number density. These uncertainty values are applied as typical throughout all cases 
and will be represented by the error bars on the solutions. 

5. Site and Measurement Characteristics 

The measurements analyzed here were performed during August and September 2006 at 
NASA/GSFC (38.99N, -76.84W). This season in the vicinity of Washington DC is frequently 
characterized by hazy, high humidity conditions and thus provides an excellent opportunity for the 
study of aerosol growth due to humidification. Such conditions are also advantageous for the 
comparison of lidar with sun photometer, because the majority of aerosols are concentrated inside the 
planetary boundary layer (PBL) thus implying that particles are reasonably spherical and that the fine 
mode dominates the PSD (Dubovik et al. 2002). 

For studying particle variation within the PBL it is desirable to perform the measurements 
starting from low altitudes. Therefore, as mentioned before, the lidar system was operated at an angle 
of 1 8 degrees with respect to the horizontal so as to minimize the vertically-resolved height to which 
the measurements were influenced by the lidar overlap function. Using this approach, backscattering 
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and extinction coefficients could be calculated starting from an altitude as low as 0.5 km above ground 
level. Measurements were performed during the nighttime, between 9-10 pm local time (01:00-02:00 
UT), which was approximately 2 hours later than the last sun photometer measurement and 
approximately 1 hour past twilight. 

Particle behavior as a function of RH is dependent upon composition and may range from 
hydrophobic to strongly hygroscopic with monotonic (smoothly varying) or deliquescent (step change) 
growth. Ambient aerosols typically found in the NASA/GSFC area consist mainly of sulfate material 
(Dubovik et al. 2002; Chen et al. 2002) and therefore significant growth of the particles can be 
expected for RH > -80%. Laboratory studies of sulfate-containing particles (Tang 1996) and field 
studies in the vicinity of NASA/GSFC bear this out (Kotchenruther et al. 1999; Fitzgerald et al. 1982). 
The value of /(RH) (defined as the ratio of total scattering at RH = 85% and 40%, respectively) for 
typical east coast regional haze aerosol has been found to be approximately 2.0 (Kotchenruther et al. 
1999; Fitzgerald et al. 1982). Increases in aerosol size over the same humidity range are on the order of 
1.3 (Fitzgerald et al. 1982). 

Profiles of atmospheric density, temperature and relative humidity (RH) were obtained from 
radiosonde measurements acquired at Dulles Airport, which is located approximately 50 km west of 
the lidar site. The radiosonde water vapor mixing ratio was calculated from these data using the 
Wexler formulation for saturation vapor pressure over a liquid surface (Wexler 1976). The radiosondes 
were launched at 00 UT, which was approximately 1-2 hours prior to the time of the lidar 
measurements. The radiosonde-derived water vapor mixing ratio data were used for the absolute 
calibration of the lidar water vapor mixing ratio measurements, which were calculated using the 
traditional Raman lidar technique (Whiteman et al. 1992). The relative humidity used in the analysis of 
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hygroscopic growth was calculated using the temperature data from the radiosonde and water vapor 
mixing ratio derived from the lidar measurements. 

5. Data selection techniques and day-by-day examples 

Multi-wavelength lidar retrievals of aerosol properties are considered here on three days with 
differing relative humidity conditions: one day with low RH (August 16, 2006) and two days with high 
RH (August 27 and September 12, 2006), and consider, for these latter two cases, the possibility that 
the measurements reveal growth of a fixed aerosol mixture due to the hygroscopic uptake of water. A 
necessary condition for studying hygroscopic growth in the atmosphere using lidar is that the same 
aerosol mixture exists at all altitudes being studied. This is achieved when the atmosphere is well- 
mixed through convection or if similar aerosols are transported to the measurement site. With similar 
aerosols being present throughout a profile, any changes in particle size due to changes in RH can be 
attributed to water uptake or release by the aerosols and hygroscopic growth curves can potentially be 
calculated. 

We used a variety of tools to detennine whether well-mixed conditions of aerosols were 
expected to exist within the lidar profile. Back trajectory analysis was used to determine the source 
region of the aerosols as a function of altitude. If the source region for the aerosols was independent of 
altitude for a given case, we took this to be evidence of similar aerosol type existing in the lidar 
profiles. We also took the presence of an altitude-independent lidar-derived water vapor mixing ratio 
to be an indication of well-mixed conditions. Finally, we investigated the derived number density from 
the lidar itself. If this did not vary beyond the uncertainty bars of the retrieval, we also took this to be 
an indication of well-mixed conditions. We considered all three of these tests for the cases to be 
presented in the next section. As we will show, two of the three datasets, those of August 27 and 
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September 12, 2006, were found to present a fixed aerosol mixture in a region of the profile. However 
only the latter case also provided the range of RHs needed to study hygroscopic growth. The weakness 
of the analysis presented here is that the temperature data used to convert the lidar mixing ratio into 
RH was provided by a radiosonde launched approximately 50 km from the lidar site. This introduces 
sufficient uncertainty in the temperature profile at the lidar site that the hygroscopic growth curves 
calculated here can only be taken as indicative of potential of the multi-wavelength technique. 

5.1. Measurements From August 16, 2006 

This day was characterized by low humidity and the absence of clouds. Fig. 3a shows the 
vertical profiles of water vapor mixing ratio obtained from the lidar and radiosonde measurements. The 
shapes of the profiles are similar, supporting the use of the sonde data for calibration of the lidar water 
vapor mixing ratio. It also suggests that the meteorological conditions were similar between the 
location of the sonde launch and the lidar site. From this point forward in the manuscript, profiles will 
be shown with respect to the vertically resolved altitude h. 

The lidar-derived RH profde agrees quite well with the sonde data, as is shown on Figure 3a. 
The RH on 16 August was below 65% in the height interval of 1.0- 1.8 km, so enhancement of the lidar 
backscatter due to particle hygroscopicity is expected to be small. The altitude dependence of aerosol 
properties measured and derived from the lidar on 16 August 2006 is shown in Fig.3b-e. The top of the 
boundary layer is at ~1.9 km and it coincides with the drop of water vapor content. Both extinction and 
backscattering have a secondary maximum near the PBL top likely due to accumulation of aerosols. 
The Angstrom exponent is found to rise with height in the PBL, which suggests that higher altitudes 
are dominated by smaller particles. This is confirmed by Fig. 3c, showing the vertical profiles of mean 
and effective radius derived from the lidar data. In the interval between 1.0 and 1.8 km, the retrieved 
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mean radius indicates a decrease in particle size. The lack of particle size increase with height is 
consistent with RH conditions below the threshold for hygroscopic growth. Fig. 3d shows the retrieved 
number and volume densities on August 16. From Fig. 3d we conclude that the previously mentioned 
maximum in particle extinction at ~1.8 km appears to be more likely due to an increase in particle 
number density than in particle size. 

Because August 16 was cloudless it provided a good opportunity to compare the lidar and sun 
photometer retrievals. The column- integrated particle parameters derived from AERONET and lidar 
measurements are summarized in Table 1. To derive column- integrated particle parameters from lidar 
data, we summed the PSDs at each height layer assuming that the PSD was constant below the lowest 
lidar retrieval. From this column-integrated PSD mean, the effective radius and volume density were 
calculated. To get column- averaged refractive indices, the values of m R and mi in each layer were 
weighted by the volume density of the layer (weighting by extinction yielded similar results): 

h 

“max 

f V(h)m(h)dh 

■=-S . 

| V(h)dh 

0 

Aerosol optical thic kn esses provided by AERONET at 355 and 532 nm are 0.32±0.015 and 
0.18±0.009. The corresponding values from lidar are 0.3±0.03 and 0. 1 8+0.01 8. The agreement is good 
keeping in mind that we extrapolated the extinction by assuming that the lowest retrieved lidar 
extinction value at 500m was constant down to the surface. The effective radius r e ff obtained from 

inversion of sun photometer measurements was r eff =0 .22±0.02 pm for the total particle size 
distribution (PSD) and r eff=0 . 12+0.01 pm for the fine mode, while the lidar retrievals for the total 
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PSD show effective radius decreasing with height from 0.28 pm to 0. 12 jam with an uncertainty of 0.04 
- 0.06 jam. The effective radius derived from the height integrated PSD is 0.21±0.04 jam for the total 
PSD and 0.15±0.03 pm when only fine mode is considered. The height profiles of real ihr and 
imaginary mi parts of refractive index are shown in Fig.3e. The uncertainty of the refractive index 
retrieval is estimated to be ±0.05 for the real part and 50% for the imaginary part. On 16 August ihr 
rises with altitude from 1.4 to 1.5, while mi stays constant at -0.01 (within the retrieval uncertainty). 
The same figure shows m R and mi obtained from the sun photometer. The spectral variation of the 
indices is small according to the sun photometer retrievals: m R changes from 1.44 to 1.43 and mi from 
0.0087 to 0.0094 in the 440-1020 nm spectral range supporting the use in this case of a spectrally 
invariant index of refraction in the lidar retrievals. Height averaged values of refractive index are 

m R = 1 45 anc * mi = 0.01 > which are close to the results from AERONET. From this comparison we 

can conclude that the height-averaged values for both effective radius and refractive index derived 
from the lidar are in reasonable agreement with the results of sun photometer measurements. 

Fig.4 shows the column-integrated PSD provided by the sun photometer. The distribution is 
bimodal with maximum at 0.1 1 pm for fine and at 5 pm for coarse mode. The same figure shows the 
lidar derived PSDs at the altitudes of 1.1 and 1.7 km. The fine modes retrieved by the two instruments 
are in reasonable agreement, but the lidar-derived PSD begins to decrease for particles slightly larger 
than 1 urn while the sun photometer PSD does not peak until -5 urn. This result could be influenced by 
the fact that the size retrievals shown in fig 4c are increasing as height decreases, implying that larger 
particles may exist below the 0.5 km lower limit of the lidar retrievals. Another influence is the fact 
that the lidar retrievals, which use a maximum wavelength of 1064 nm (-1 um), have decreasing 
sensitivity to particles with sizes greater than -2 um. Simulations (Veselovskii et al. 2004) have shown, 
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however, that the estimation of coarse mode parameters is possible although additional information 
such as would be provided through the use of longer wavelengths would be beneficial to the retrieval. 

The columnar volume density from sun photometer was 0.056±0.005 pm /pm“ for total PSD 

3 2 

and 0.029±0.004 pm /pm for the fine mode. The corresponding lidar-derived values are 0.035 ±0.007 
pm /pm" and 0.025 ±0.005 pm /pm . In accordance with Fig.4 the main part of the lidar-derived 
density was contained in the fine mode and the corresponding value agrees well with the sun 
photometer measurements. 

5.2. Measurements from 27 August 2006 

The measurements on August 27, 2006 were performed under RH conditions sufficiently high 
to potentially support hygroscopic particle growth studies. At the time of the measurements, however, 
the back trajectories suggest that three separate air masses were influencing the vertical profile of 
aerosols present at the lidar site: 1) for the height region of 800-2000 m, the origin of the air was 
generally southerly, 2) for the height region of 2350-2800 m, the origin was generally south-westerly 
from the Virgina -Kentucky region and 3) for the height region of 2950-3400 m, the origin is also to the 
southwest but from the Alabama-Georgia region. Therefore, without sufficient mixing we would not 
expect to find the same type of aerosols at all altitudes thus preventing a clean determination of 
hygroscopic growth characteristics even if co-located radiosonde temperature data had been available. 
It should be noted that clouds were present on this day, which prevented any sun photometer retrievals 
from occurring so no comparisons with lidar were possible. 

Fig. 5 shows the retrievals for August 27. The top of the boundary layer is at ~3.5 km, the 
radiosonde and lidar measured RH profiles are quite similar and show relative humidity increasing with 
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altitude from values below 50% at the bottom of the lidar retrieval around 500m to above 90% for 
altitudes between 2.5 - 3.5 km. The water vapor mixing ratio profile shows a slight decrease with 
altitude suggesting that the boundary layer may have been well-mixed during the daytime but that, at 
the times and locations of both the sonde launch and the lidar measurements, mixing was not complete 
in the boundary layer. The incomplete mixing present this day does not make this a good case for the 
study of hygroscopic growth. Nonetheless, the retrievals reveal a strong correlation between particle 
properties and RH. Inside the altitude range 0.5- 1.8 km the particle size decreased slightly with height 
(Fig 5c). At an altitude of 1.8 km, RH exceeded 80% and the Angstrom exponent at this height starts to 
decrease from 1.5 to 0, which indicates an increase in particle size as shown in Fig. 5c. A decrease in 
the Angstrom exponent is accompanied by an increase in the particle extinction and backscattering 
coefficients. The mean and effective radii grow rapidly above 2.0 km reaching 0.30 pm and 0.6 pm 
correspondingly. The fast growth of effective radius and volume density in the altitude interval 2.0 - 
3.3 km indicates the appearance of big particles resulting in an increase of extinction. The retrieved hir 
on August 27 was very close to that of pure water and varied from 1.32 to 1.36 in the height range 
considered. This variation in hir can be considered to be within retrieval uncertainty. The imaginary 
part decreases above 2.0 km as shown in Fig.5e. This could be attributable to water uptake in the 
regions of the aerosol population characterized by higher RH. 

5.3. Measurement from 12 September 2006 

In our measurement sessions, high humidity was usually observed during cloudy periods (e.g., 
August 27 (described above), so there were few opportunities to make a comparison between the lidar 
retrievals and those of sun photometer. One case where this was possible, however, occurred on 
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September 12, 2006. The profiles of lidar-derived particle parameters for this case are shown in Fig. 6 
and the comparison of lidar and AERONET derived particles parameters is summarized in Table 2. 

The back trajectory analysis for September 12 suggests that for the altitude regions of 1150 - 
2500 m at the lidar site, the origin of the aerosols was from a common region to the NW of the lidar 
site. We take this to imply that the aerosol mixture between these two altitudes is likely to be the same 
at the lidar site. This coupled with the wide range of RH present makes this case, of the three presented 
here, the one most likely to show hygroscopic growth behavior similar to more controlled in-situ 
studies and thus pennit hygroscopic growth parameters to be calculated. 

The RH shown in Fig 6 on September 12 exhibited an increase with altitude reaching a 
maximum at 1.5 km and then decreasing beyond 1.8 km. The radiosonde data indicate saturated 
conditions between the altitudes of 1.2 and 2.0 km revealing the likely presence of a cloud. The lidar 
data, however, did not indicate the presence of a cloud illustrating again the desire to have co-located 
radiosonde data for studies such as these. 

The altitude interval between 1.2 and 2.0 km is characterized by an increase in particle 
extinction from 0.05 to 0.4 km' 1 at 355 nm and by a decrease in the Angstrom exponent from 1.2 to 0.5 
(Fig. 6b). The Angstrom exponent derived from the sun photometer at 340 nm and 500 nm wavelengths 
is approximately 0.9, while lidar derived value varies from 0.4 to 1.8. This is another example 
illustrating that significant changes in particle parameters can occur as a function of altitude that are 
masked by a column-integrated retrieval from the sun photometer. 

The effective radius varies from 0.13 to 0.5 pm while the column- integrated value from the sun 
photometer is 0.3±0.02 pm. The columnar volume concentrations from the sun photometer when the 
total PSD and only fine mode is considered are 0.073±0.007 pm 3 /pm 2 and 0.05±0.0075 pm 3 /pm 2 
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3 2 

respectively. The corresponding lidar-derived values are of 0.076±0.015 pm /pm and 0.057±0.012 

3 2 

pm /pm“. Thus the results from both instruments are in reasonable agreement. 

The real part of lidar-derived refractive index does not change much with altitude and the 

altitude-averaged value ntR -1.34±0.05 agrees well with the sun photometer result 1.34±0.04. The 

imaginary part mi derived from lidar varies with height between 0.0 1 5+0. 0075 and 0.007±0.0035. The 

height averaged value of imaginary part derived from lidar is mi = 0.009 ±0.0045 > while for the sun 

photometer mi=0.05±0.0075 > thus these values are in agreement. 

We also studied the correlation between the particle size and the Angstrom exponent. Fig. 7 
shows the dependence of effective radius on Angstrom exponent for the ensemble of data. There is an 
inverse dependence of k on r e ff that is nearly linear, implying that the lidar-derived k can be used for a 
preliminary estimation of particle size. The same figure shows the results for AERONET 
measurements on 16 August and 12 September. These results are in good agreement with lidar 
measurements. 

6. Potential of multi-wavelength lidar for aerosol hygroscopic growth studies. 

As discussed in the earlier sections, the lidar measurements on September 12 were acquired 
under conditions considered to be supportive of hygroscopic growth (range of RH and similar aerosol 
type in the profile) while the data from August 27 were not. We will now consider the impact of RH on 
the particle parameters for these two cases taking the radiosonde temperature data to represent the 
temperature structure at the lidar site. The uncertainty in this assumption, due to the separation of the 
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lidar and radiosonde sites, implies that the hygroscopic growth curves calculated from the lidar data 
should be taken as illustrating the potential of the technique and not as absolute quantifications. 

The change of aerosol scattering properties as a function of RH is usually described by the 
humidification factor (or hygroscopic growth factor) /(RH) used to characterize the response of a 
homogeneous mixture of aerosols that is subjected to changing values of RH. It is defined as the ratio 
between wet and dry scattering cross sections as follows: 


RRH) = 


a(RH ) 
a(RH ref ) 


, where a(RH re f) corresponds to the extinction coefficient at a low, reference RH 


value, which is obtained at some altitude level. The humidification factor /RH) can be parameterized 
using the Hanel model (Hanel 1976), as it was done by many investigators (Kotchenruther et al. 1998, 
1999; Raut and Chazette 2007; Randriamiarisoa et al. 2006). 

a = a lby (l-RH)~ g (13) 


where oidry is the extinction of dry particles. For the reference value of RH this relationship can be 
rewritten as: 



(14) 


We now calculate f(RH) from the data of August 27, when different source regions were 
identified in the data and mixing was incomplete, and September 12 where the measurement conditions 
are considered supportive of hygroscopic growth studies. The results for August 27 are shown in Fig 8 
where a reference RH value of 60% was used. A general increase in/ is observed with increasing RH, 
however the curve shows a great deal of scatter and differs significantly from the Hanel model results 
for g=0.9, which is the best fit for 12 September results. We take the strong departure of the 
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experimental curve from the simulated curve shown in fig 8 to reflect the fact that the conditions 
present on August 27 were not supportive of hygroscopic growth studies: different source regions of 
the aerosols were identified along with incomplete mixing indicating the likely presence of different 
aerosol types as a function of altitude in the profile. However, the general trend of increasing extinction 
as RH increases is consistent with the retrievals shown in figures 5 indicating that changing RH 
strongly influenced the profile of retrieved parameters. 

Considering now data from September 12, when conditions were more suitable for hygroscopic 
growth studies, Fig. 9 shows the dependence of extinction coefficient at 532 nm on RH for three 
altitude ranges. The humidification factor /(RH) is also shown in figure 9. Again, the reference value 
RH re f was chosen to be 60%, because lower values of RH were not accessible in these measurements. 
The value of/fRH) at RH=80% in our measurements is 1.78 at 532 nm. Using Figure 13 in Fitzgerald 
et al. 1982, an /( R H ) at A,— 550 nm for the same humidity range is —1 .45. Using Figure 5 in 
Kotchenruther et al. 1 999, /(RH) values of ~1.5 and -1.7 are obtained for clean and anthropogenic- 
influenced air respectively over the same RH range. Thus the obtained values from the lidar are 
reasonably close to previous results for typical regional haze aerosol on the US East coast. 

The best fit to experimental data at 532 nm for RH<85% is achieved at g=0.9, the 
corresponding curve for which is shown in Fig. 9. The derived value of g is in reasonable agreement 
with earlier studies (Randriamiarisoa et al. 2006) showing variation of this parameter between 0.27 and 
1.35, so the derived value of g looks reasonable. For RH exceeding 85% the lidar-derived 
humidification factor rises faster than the model prediction. This may be due to deliquescent growth 
implying that the simplified parameterization used in expression (14) is insufficient to describe the 
physics involved and that more complicated models may be needed (Kotchenruther et al. 1999). The 


23 



hysteresis behavior displayed in Fig 9 may also indicate the presence of deliquescent compounds in the 
aerosols (Randriamiarisoa et al. 2006). 

The multi-wavelength lidar technique also allows estimation of aerosol size growth factor 


f r (RH) = 


r(RH) 

r{RH,. ef ) 


, which is the ratio between the wet and the quasi-dry particle radius. The 


dependence of the aerosol mean radius on RH together with size growth factor is shown in Fig. 10. For 
comparison the same picture shows the hygroscopic growth of ammonium sulfate (NFL^SC^, 
calculated for the dry particle radius 0.1 pm using the expressions from (Tang 1996). Although the 
scattering of data in Fig. 10 is significant, the value of f( 85%) is estimated to be about 1.3 which is a 
typical number for the east coast regional haze aerosol (Kotchenruther et al. 1999). 

The retrieval of integral particle parameters, such as volume density, is more stable than 
number density retrievals. Therefore it is interesting to consider the volume growth factor f,(RH) that 
may be more appropriate for the study of the hygroscopic growth process using multi-wavelength 


Raman lidar: X 


V(RH) 

V(RH ) ' defined this way, f v (RH) represents the amount of water accepted in 


the process of particle growth. The corresponding results are shown in Fig. 1 1 . This factor reveals a 
smoothly varying dependence on RH, which is due to larger range of volume variation to compare with 
radius. For comparison, the results from August 27 are also shown. For that day V(RH) does not follow 
the model curve which could be another indication of the presence of different aerosol types along the 
sounding path. 


7. Summary and conclusions 
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The main intention of this paper was to demonstrate the ability of the multi-wavelength Raman 
lidar technique to profile the vertical distribution within the planetary boundary layer of critical aerosol 
parameters such as mean and effective radius, Angstrom exponent, complex refractive index, number 
and volume densities under varying relative humidity conditions. A newly designed multi-wavelength 
Raman lidar and upgraded inversion algorithm pennits the simultaneous retrieval of microphysical 
particle parameters at all height in the profile. Regardless of the degree of mixing, the source of 
aerosols or other considerations, we found a strong positive correlation between RH and aerosol 
particle extinction, size and volume. To support the results of the lidar measurements, comparisons 
with column-averaged quantities from a co-located AERONET sun photometer data were performed. 
The lidar measurements of fine mode parameters agreed well with the sun photometer, but at present, 
the lidar technique has difficulty in fully representing the coarse mode. A contributing factor could be 
that the particles exceeding 2 pm may have been concentrated at altitudes lower than were retrieved by 
lidar in this study (0.5 km), hence our next goal is to decrease the minimum useful altitude of the 
measurements. These comparisons make particularly clear the need to consider the possibility of 
hygroscopic growth when interpreting AERONET column-averaged data. 

Three datasets were studied in detail. One was characterized by low RH conditions where no 
evidence of particle growth with increasing RH was observed and two were acquired under higher RH 
conditions where a strong correlation was observed between aerosol size and RH. Considering these 
two latter cases, hygroscopic growth factors were calculated. However the first of the two cases on 
August 27 was likely characterized by differing aerosol mixtures in the vertical while the second case, 
that on September 12, 2006, was judged to have both a range of RH and a fixed aerosol mixture in a 
portion of the profile making it suitable for calculations of hygroscopic growth parameters. The 
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calculated hygroscopic growth curves for these two cases differed considerably with the latter case on 
September 12 being more consistent with published in literature. 

The results presented here demonstrate the potential of the multi-wavelength lidar technique for 
the study of aerosol hygroscopic growth in the atmosphere. Although the derived hygroscopic growth 
parameters agree well with modeling, they must be considered preliminary due to the lack of 
temperature data at the lidar site and the inability of fully resolve the coarse mode of the PSD. For 
validation of the MW technique, the temporal and spatial separation between lidar and radiosonde 
measurements should be minimized and cases should be selected where the PSD is dominated by the 
fine mode. 
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Tables 


Table 1. Column integrated aerosol parameters derived from lidar and AERONET on 16 August 


2006. 


Parameter 

AERONET 

Lidar 

AOT 

355 nm 

0.32±0.015 

0.3±0.03 

532 nm 

0.18±0.009 

0. 1 8+0.01 8 

V tot (pm 3 / pm 2 ) 

0.056±0.005 

0.035±0.007 

V fine (pm 3 / pm 2 ) 

0.029±0.004 

0.025±0.005 

tot , N 

r eff (pm) 

0.22±0.06 

0.21±0.04 

4 ne (pm) 

0.123±0.035 

0.15±0.03 

m R 

1.44±0.04 

1.45±0.05 

mi 

0.01±0.003 

0.01±0.005 


3 




Table 2. Column integrated aerosol parameters derived from lidar and AERONET on 
12 September 2006. 


Parameter 

AERONET 

Lidar 

AOT 

355 mn 

0.38±0.019 

0.41±0.08 

532 nm 

0.27±0.013 

0.31±0.06 

V tot (pm 3 /pm 2 ) 

0.073±0.007 

0.076±0.015 

V fme (pm 3 /pm 2 ) 

0.05±0.0075 

0.057±0.012 

tot , V 

r eff (pm) 

0.301±0.09 

0.33±0.07 

r% e (pm) 

0.21±0.06 

0.25±0.05 

m R 

1.37±0.04 

1.34±0.05 

mi 

0.005±0.0015 

0.009±0.0045 
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Captures to the figures 

Fig.l. Lidar signal at 1064 nm measured on 2 September, 2006 together with simulated Raleigh 
signal (dash-dot). 

Fig. 2. Vertical profiles of aerosol backscattering (solid) and extinction (dash-dot) coefficients at 
355, 532 and 1064 nm measured on 22 August 2006. 

Fig 3. Vertical profiles of aerosol parameters measured on 16 August: (a) water vapor mixing ratio 
Wm and relative humidity RH derived from the lidar data (solid lines) and from the sonde 
measurements (line+symbol); (b) aerosol extinction coefficients at 355 nm together with Angstrom 
exponent k\ (c) mean r mean and effective r eff radius, dashed line shows the effective radius for the 

fine mode ( r e ff e ) and for the total size distribution (kff ) obtained by AERONET; (d) number N and 
volume V density; (e) real m R and imaginary mi part of refractive index, dotted lines show the 
results from AERONET. 

Fig. 4. Particle size distribution derived from lidar measurements at altitudes 1.1 and 1.7.km on 16 
August. Open symbols show the columnar PSD from AERONET. 

Fig. 5. Vertical profiles of aerosol parameters measured on 27 August: (a) water vapor mixing ratio 
Wm and relative humidity RH derived from the lidar data (solid lines) and from the sonde 
measurements (line+symbol); (b) aerosol extinction coefficients at 355 nm together with Angstrom 
exponent k\ (c) mean r mean and effective r e ff radius, dashed line shows the effective radius for the 

fine mode (kff 6 ) and for the total size distribution (kff ) obtained by AERONET; (d) number N and 
volume V density; (e) real m R and imaginary mi part of refractive index, dotted lines show the 
results from AERONET. 

Fig. 6. Vertical profiles of aerosol parameters measured on 12 September : (a) water vapor mixing 
ratio Wm and relative humidity RH derived from the lidar data; (b) aerosol extinction coefficients at 
355 nm together with Angstrom exponent k: (c) mean r mean and effective r e ff radius, dashed line 
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shows the effective radius for the fine mode (r eff ) and for the total size distribution ( r cff ) obtained 
by AERONET; (d) number N and volume V density; (e) real m R and imaginary mi part of refractive 
index, dotted lines show the results from AERONET. 

Fig. 7 Correlation between lidar derived Angstrom exponent and effective radius. Solid stars 
represent the results for AERONET on 16 August and 12 September. 

Fig. 8. /(RH) from the dataset of August 27, 2006 for the ranges 800 - 2000 m (stars) and 2000 - 
2600 m (squares). For these ranges two distinct source regions for aerosols were identified by back 
trajectory analysis. Solid line represents the parameterization fit of/(RH) for g=0.9. 

Fig. 9. Dependence of extinction coefficient at 532 nm on RH and humidification factor /(RH) 
derived from 12 September data. Results are shown for altitude ranges: 1.0-1. 3 km (solid squares), 
1.3 - 1.8 km (open stars) and 1.8-2. 5 km (open circles). Solid line represents the parameterization fit 
of^(RH) for g=0.9. 

Fig. 10. Dependence of aerosol mean radius on RH derived from 12 September 2006 data together 
with size growth factor A(RH). Solid line shows the hygroscopic growth of (NH 4 ) 2 S0 4 particle with 
r 0 =0.1 pm. 

Fig. 11. Dependence of aerosol volume concentration on RH derived from 12 September 2006 data 
together with volume growth factor. For comparison the results from 27 August, when the condition 
of well mixing are not fulfilled are also shown. 
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Fig. 1 . Lidar signal at 1064 nm measured on 2 September, 2006 together with simulated Raleigh 
signal (dash-dot). 
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Fig. 2. Vertical profiles of aerosol backscattering (solid) and extinction (dash-dot) coefficients at 
355, 532 and 1064 nm measured on 22 August 2006. 
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Fig. 3. Vertical profiles of aerosol parameters measured on 16 August: (a) water vapor mixing ratio 
Wm and relative humidity RH derived from the lidar data (solid lines) and from the sonde 
measurements (line+symbol); (b) aerosol extinction coefficients at 355 mn together with Angstrom 
exponent k\ (c) mean r mean and effective r e ff radius, dashed line shows the effective radius for the 

fine mode ( r e ff e ) and for the total size distribution ( r ) obtained by AERONET; (d) number N and 
volume V density; (e) real mR and imaginary mi part of refractive index, dotted lines show the 
results from AERONET. 
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Fig. 4. Particle size distribution derived from lidar measurements at altitudes 1.1 and 1.7. km on 16 
August. Open symbols show the columnar PSD from AERONET. 


4 




Altitude, km Altitude, km 


Vapor mixing ratio, g/kg 



RH, % 


Angstrom exponent 


0.0 0.5 1.0 1.5 2.0 



Extinction, km ' 1 


4 




Altitude, m Altitude, km 



V (fim 3 /cm 3 ) 

o 20 40 60 80 

3.5 





m l 

0.000 0.005 0.010 0.015 0.020 



Fig. 5. Vertical profiles of aerosol parameters measured on 27 August: (a) water vapor 
mixing ratio Wm and relative humidity RH derived from the lidar data (solid lines) and 
from the sonde measurements (line+symbol); (b) aerosol extinction coefficients at 355 
mn together with Angstrom exponent k: (c) mean r mean and effective r e ff radius, dashed 

line shows the effective radius for the fine mode ( faff”) and for the total size distribution 

( hff ) obtained by AERONET ; (d) number N and volume V density; (e) real m R and 
imaginary mi part of refractive index, dotted lines show the results from AERONET. 
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Fig. 6. Vertical profiles of aerosol parameters measured on 12 September : (a) water vapor mixing ratio 
Wm and relative humidity RH derived from the lidar data; (b) aerosol extinction coefficients at 355 mn 
together with Angstrom exponent k; (c) mean r mean and effective r e ff radius, dashed line shows the 

effective radius for the fine mode ( r e f ™ e ) and for the total size distribution ( r^ ) obtained by 
AERONET; (d) number N and volume V density; (e) real mR and imaginary mi part of refractive 
index, dotted lines show the results from AERONET. 
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Fig. 7. Correlation between lidar derived Angstrom exponent and effective radius. Solid stars represent 
the results for AERONET on 16 August and 12 September. 
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Fig. 8.y(RH) from the dataset of August 27, 2006 for the ranges 800 - 2000 m (stars) and 2000 - 2600 
m (squares). For these ranges two distinct source regions for aerosols were identified by back trajectory 
analysis. Solid line represents the parameterization fit of/(RH) for g=0.9. 
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Fig 9. Dependence of extinction coefficient at 532 nm on RH and humidification factor /(RH) derived 
from 12 September data. Results are shown for altitude ranges: 1.0- 1.3 km (solid squares), 1.3 - 1.8 km 
(open stars) and 1.8-2. 5 km (open circles). Solid line represents the parameterization fit of /(RH) for 
g=0.9. 
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Fig. 10. Dependence of aerosol mean radius on RH derived from 12 September 2006 data together with 
size growth factor /i(RH). Solid line shows the hygroscopic growth of (NFL^SC^ particle with r 0 =0.1 
pm. 
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Fig. 11. Dependence of aerosol volume concentration on RH derived from 12 September 2006 data 
together with volume growth factor. For comparison the results from 27 August, when the condition of 
well mixing are not fulfilled are also shown. 
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